Excitation spectra, spin structures, and entanglement characteristics 
of four-electron double-quantum-dot artificial molecules 



o 
o 

(N 

in 



Ying Li, Constantine YannouleasQ and Uzi LandmarQ 
School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332-0430 

(Dated: 30 January 2009) 

Energy spectra, spin configurations, and entanglement characteristics of a system of four electrons 
in lateral double quantum dots are investigated using exact diagonalization (EXD), as a function 
of interdot separation, applied magnetic field, and strength of interelectron repulsion. A distinctly 
different quantum behavior is found compared to that of circular single quantum dots. As a function 
of the magnetic field, the energy spectra exhibit a low-energy band consisting of a group of six states, 
with the number six being a consequence of the conservation of the total spin of the four electrons 
and the ensuing spin degeneracies. These six states appear to cross at a single value of the magnetic 
field, with the crossing point becoming sharper for larger interdot distances. As the strength of the 
Coulomb repulsion increases, the six states tend to become degenerate and a well defined energy 
gap separates them from the higher-in-energy excited states. 

The appearance of the low-energy band is a consequence of the formation of a Wigner super- 
molecule, with the four electrons (two in each dot) being localized at the vertices of a rectangular 
parallelogram. Using the spin-resolved pair-correlation functions, one can map the EXD many-body 
wave functions onto the spin functions associated with the four localized electrons. Analogies with 
nanomagnets, such as finite Heisenberg clusters, are discussed. The ability to determine associated 
spin functions enables investigations concerning entanglement properties of the system of four elec- 
trons. In particular, the formation of Wigner supermolecules generates EXD solutions belonging to 
the class of strongly entangled states referred to as A r -qubit Dicke states in the quantum-information 
literature. 

PACS numbers: 73.21. La, 31.15.V-, 03.67.Mn, 03.65.Ud 
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I. INTRODUCTION 

The field of two-dimensional (2D) semiconductor 
quantum dots (QDs) has witnessed rapid expansion 
in the last several years, both experimentally^ and 
theoretically.—^ Along with fundamental interest in the 
properties of such systems, and as a test ground for highly 
correlated electrons, a major motivation for these grow- 
ing endeavors has been the promising outlook and po- 
tential of quantum dots concerning the implementation 
of solid-state quantum computing and quantum informa- 
tion devices£i£i£ To this effect highly precise control of 
the space and spin degrees of freedom of a small num- 
ber N of confined electrons (down to an empty^ ' 10 : 11 
QD) needs to be achieved, and experimentally this was 
demonstrated recently for two electrons in a lateral dou- 
ble quantum dot molecule (see Ref. 0, and references 
therein). From the theoretical standpoint, high-level 
computational methods that reach beyond the level of 
mean-field approximation are needed; 5 - with the ability 
to provide solutions that preserve all the symmetries of 
the many-body Hamiltonian. For example, conserva- 
tion of total spin is essential for describing spin entan- 
glement in small assemblies of electrons that exhibit lo- 
calization in space. Such assemblies may be viewed as 
finite Heisenberg spin cluster s 12 ' 13 whose quantum be- 
havior (due to finite-size fluctuations and correlation ef- 
fects) differs drastically from the behavior expected from 
magnetic systems in the thermodynamic Unfit i 12 ' 

There is an abundance of experimental and theoreti- 
cal publications concerning circular single quantum dots 



with a small number of electrons ! 1 ' 3 ' 4 ' 5 ' 15 ' 16 ' 17 In this pa- 
per, we use exact diagonalization 5 (EXD) to investigate 
the properties of lateral double quantum dots (DQDs) 
containing four electrons. DQDs are referred to also as 
artificial molecules. DQDs containing two electrons have 
been already studied extensively both experimentally 2 
and theoretically, 5 . ! 18 ! 19 However, experimental studies of 
DQDs with three or four electrons have not been per- 
formed as yet; we are aware of a single theoretical study 
of a lateral DQD with three electrons^ and another 
one of two laterally coupled quantum rings with three 
electrons <2I 

In light of the novel quantum behavior discovered in 
our investigations (compared to circular QDs, both con- 
cerning the spectra and entanglement aspects), we hope 
that the present work would serve as an impetus for fur- 
ther experimental studies on lateral DQDs. In particular, 
as a function of the magnetic field, we find that: (1) A 
low-energy band of six states develops as the strength 
of the Coulomb repulsion increases, separated by an en- 
ergy gap from the other excited states, and (2) All six 
states appear to "cross" at a single value of the magnetic 
field. The crossing point gets sharper for larger interdot 
distances. We find that the specific number of crossing 
states (six) derives from the spin degeneracies and mul- 
tiplicities expressed by the branching diagram^ 

The formation of the low-energy band is a consequence 
of the localization of the four electrons within each dot 
(with two electrons in each dot). This localization leads 
to formation (with increasing strength of the Coulomb 
repulsion) of a Wigner supermoleculej^ 3 . with the four lo- 
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calized electrons (two in each dot) being located at the 
corners of a rectangular parallelogram. Using the spin- 
resolved pair-correlation functions, we show how to map 
the EXD many-body wave functions onto the spin func- 
tions associated with the four localized electrons. This 
mapping leads us naturally to study analogies with finite 
systems described by a model Heisenberg Hamiltonian 
(often referred to as finite Heisenberg clusters). We fur- 
ther discuss that the determination of the equivalent spin 
functions enables investigations concerning the entangle- 
ment properties of the EXD solutions. In particular, we 
show that the formation of Wigner supermolecules leads 
to strongly entangled states known in the literature of 
quantum information as TV-qubit Dicke states i 24 ' 25 ^ 26 ' 27 
We finally mention that the trends in the excitation 
spectra (e.g., formation of a low-energy band) and en- 
tanglement properties (e.g., mapping to spin functions 
of localized electrons) found in the case of a double 
quantum dots have many analogies with those found in 
other deformed configurations, and in particular single 
anisotropic quantum dots; see, e.g., the case of three elec- 
trons in Ref. Hfl. 



II. TWO-DIMENSIONAL 
TWO-CENTER-OSCILLATOR CONFINING 
POTENTIAL 

In the two-dimensional two-center-oscillator (TCO), 
the single-particle levels associated with the confining po- 
tential of the artificial molecule are determined by the 
single-particle hamiltonian 2 ^ 

H = T + l -m*uy + \m*Ji h x'l 

+ V neck (x) + h k + ^-B-s , (1) 
n 

where x' k — x — x k with k = 1 for x < (left) and 
k = 2 for x > (right), and the h k 's control the relative 
well-depth, thus allowing studies of hetero-QDMs. y de- 
notes the coordinate perpendicular to the interdot axis 
(x). T = (p - eA/c) 2 /2m*, with A = 0.5(-By,Bx,0), 
and the last term in Eq. (fTJ) is the Zeeman interaction 
with g* being the effective g factor, [Ib the Bohr magne- 
ton, and s the spin of an individual electron. The most 
general shapes described by H are two semiellipses con- 
nected by a smooth neck [V r „ ec fe(x)]. x± < and X2 > 
are the centers of these semiellipses, d = x 2 — x\ is the 
interdot distance, and m* is the effective electron mass. 

For the smooth neck, we use V neck (x) = 
\m*Lu 2 xk [C k x'* + V k x' k 4 }8(\x\ - \x k \), where 6{u) = 
for u > and 8(u) — 1 for u < 0. The four constants 
C k and T> k can be expressed via two parameters, as 
follows: C k = (2 - 4e\)/x k and V k = (1 - 3e b k )/x 2 , 
where the barrier-control parameters e b k = (V& — h k )/Vo k 
are related to the actual (controlable) height of the 
bare interdot barrier (VJ,) between the two QDs, and 
V ok = m*cj 2 xk x 2 /2 (for hi = h 2 , Vol = V 02 = V ). 



The single-particle levels of H, including an external 
perpendicular magnetic field B, are obtained by numer- 
ical diagonalization in a (variable-with-separation) basis 
consisting of the eigenstates of the auxiliary (zero-field) 
hamiltonian: 

p 2 1 1 

Ho = 2^ + 2 y y + 2 m * xkX * + hk • (2) 

The eigenvalue problem associated with the auxiliary 
hamiltonian [Eq. @] is separable in x and y, i.e., the 
wave functions are written as 

(fi(x,y) =X^,(x)Y n (y), (3) 

with i = {/i, n}, i — 1,2, .... K . 

The Y n (y) are the eigenfunctions of a one-dimensional 
oscillator, and the X^(x < 0) or X^(x > 0) can be 
expressed through the parabolic cylinder function a 30 ! 31 
[/[7fc, (-l) fe £ fc ], where ^ = x' kV / 2m*uj xk /h, ^ k = (-E x + 
h k )/ (hui xk ), and E x = (/i + 0.5)fuj x i + hi denotes the x- 
eigenvalues. The matching conditions at x = for the 
left and right domains yield the x-eigenvalues and the 
eigenfunctions X„(x). The n indices are integer. The 
number of fi indices is finite; however, they are in gen- 
eral real numbers. 

In the Appendix, we discuss briefly the energy spec- 
tra associated with the single-particle states of the two- 
center oscillator Hamiltonian given by Eq. |T|) . We follow 
there the notation presented first in Ref. [32. For further 
details, see Ref. l33l 

In this paper, we will limit ourselves to QDMs with 
X2 = —X\ and haj y = Twj x \ = hu X 2 — hujQ. How- 
ever, in several instances we will compare with the case 
of a single elliptic QD where x 2 — —x\ = and 
hujy hcu x = huj x i = hu X 2. In all cases, we will use 
huiQ = 5.1 meV, m* — 0.070m e (this effective-mass value 
corresponds to GaAs), and K = 50 (which guarantees 
numerical convergence 3 ^) . 

III. THE MANY-BODY HAMILTONIAN AND 
THE EXACT DIAGONALIZATION METHOD 

The many-body hamiltonian Ti. for a dimeric QDM 
comprising N electrons can be expressed as a sum of 
the single-particle part H(i) defined in Eq. |T]) and the 
two-particle interelectron Coulomb repulsion, 

N N N 2 

tt = £^)+£E^ (4) 

I— 1 I— 1 J>% J 

where k is the dielectric constant and denotes the 
relative distance between the i and j electrons. 

As we mentioned in the introduction, we will use the 
method of exact diagonalization for determining^ 3 - the so- 
lution of the many-body problem specified by the hamil- 
tonian Q. 

In the EXD method, one writes the many-body wave 
function $^ XD (ri, r 2 , ■ ■ ■ , rjv) as a linear superposition 
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of Slater determinants W (ri, r2, . . . , r/v) that span the 
many-body Hilbert space and are constructed out of the 
single-particle spin-orbitals 



Xj(x,y) = Vj(x,y)a, if 1 < j < K , 



(5) 



and 



Xj{x,y) =<Pj- K (x,y)P, if K<j<2K, (6) 
where ct{j3) denote up (down) spins. Namely 



where the operator U7y interchanges the spins of electrons 
i and j provided that their spins are different; N a and Np 
denote the number of spin-up and spin-down electrons, 
respectively. 

Of great help in reducing the size of the matrices to 
be diagonalized is the fact that the parity (with respect 
to the origin) of the EXD many-body wave function is 
a good quantum number for all values of the magnetic 
field when h\ = h%. Specifically, the rry-parity operator 
associated with reflections about the origin of the axes is 
defined as 



<T> EXD /V 



,rjv) 



C?*?(n, • • • , r N ), (7) ^ (ril r 2 , r„ r 4 ) = *^(-n, -r 2 , -r 3 , -r 4 ) 



i 



where 



nsn 



Xji(rjv) ... Xon( v n) 



(8) 



and the master index / counts the number of arrange- 
ments {ji,j2, ■ ■ ■ ,Jn} under the restriction that 1 < ji < 
]2 < ■ ■ ■ < on < 2K. Of course, q — 1, 2, . . . counts the 
excitation spectrum, with q = 1 corresponding to the 
ground state. 

The exact diagonalization of the many-body 
Schrodinger equation 



EXD 

N,q 



pEXD^EXD 



(9) 



transforms into a matrix diagonalizatiom problem, which 
yields the coefficients Cj and the EXD eigenenergies 
i?^ D . Because the resulting matrix is sparse, we im- 
plement its numerical diagonalization employing the well 
known ARPACK solver^ 

The matrix elements (ty I N \7-C\ i f>ff) between the basis 
determinants [see Eq. ||5J)] are calculated using the Slater 
rules. 36 Naturally, an important ingredient in this respect 
are the two-body matrix elements of the Coulomb inter- 
action, 

e 2 roo ,00 l 

— / / dr 1 dr 2 ^*(r 1 W*(r 2 ) 1 -Ap k (r x )yi{r 2 ), 

K J-ooJ-oo 3 |ri -r 2 | 

. (10) 

in the basis formed out of the single-particle spatial or- 
bitals </?i(r), i = 1, 2, . . . , K. In our approach, these ma- 
trix elements are determined numerically. 

The Slater determinants [see Eq. (jHJ)] conserve the 
third projection S z , but not the square S 2 of the total 
spin. However, because S 2 commutes with the many- 
body hamiltonian, the EXD solutions are automatically 
eigenstates of S 2 with eigenvalues S(S + 1). After the 
diagonalization, these eigenvalues are determined by ap- 
plying S 2 onto Q^g and using the relation 



S 2 *f 



N f3 ) 2 /i- 



(ii) 



(12) 

and has eigenvalues ±1. 

One can also consider partial parity operators V x and 
V y associated solely with reflections about the x and y 
axis, respectively; of course V xy = V X V V . We note that 

unlike Vxy, the partial parities Vx and V y are conserved 
only for zero magnetic fields (B = 0). With the two- 
center oscillator cartesian basis that we use [see Eq. fl3"))], 
it is easy to calculate the parity eigenvalues for the Slater 
determinants, Eq. ([8]). that span the many-body Hilbert 
space. Because X fJi (x) and Y n (y) conserve the partial V x 
and V y parities, respectively, one finds: 



V *i 



N 



N 



(13) 



where rrii and rij count the number of single-particle 
states associated with the bare two-center oscillator [see 
the auxiliary hamiltonian Hq in Eq. (|2[)] along the x axis 
and the simple oscillator along the y direction (with the 
assumption that the lowest states have m = and n = 0, 
since they are even states). We note again that the in- 
dex /i in Eq. ([3]) is not an integer in general, while m 
here is indeed an integer (since it counts the number of 
single-particle states along the x direction). 



IV. MANY-BODY SPIN EIGENFUNCTIONS 

For completeness and for the reader's convenience, we 
outline in this section several well established (but of- 
ten not well known) properties of the many-body spin 
eigenfunctions which are useful for analyzing the trends 
and behavior of the spin multiplicities exhibited by the 
EXD wave functions for N — 4 electrons. We stress 
here that the ability to describe spin multiplicities is an 
advantage of the EXD method compared to the more 
familiar spin-density functional approaches whose single- 
determinantal wave functions preserve only the third pro- 
jection S z of the total spin, and thus are subject to "spin 
contamination" errors. As we will discuss below, the spin 
multiplicities of the EXD wave functions lead naturally 
to formation of highly entangled Dicke states . 24 i 25 i 26 i 27 

A basic property of spin eigenfunctions is that they ex- 
hibit degeneracies for N > 2, i.e., there may be more than 
one linearly independent (and orthogonal) spin functions 
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FIG. 1: The branching diagram for the spin degeneracies. 
The total-spin quantum number S is given on the vertical 
axis, and the number of particles, N, on the horizontal one. 
The numbers inside the circles give the number, g(N,S), of 
linear independent (and orthogonal) spin functions for the 
corresponding values of N and S. 

that are simultaneous eigenstates of both S 2 and S z . 
These degeneracies are usually visualized by means of the 
branching diagram 22 displayed in Fig. [TJ The axes in this 
plot describe the number N of fermions (horizontal axis) 
and the quantum number S of the total spin (vertical 
axis). At each point (N, S), a circle is drawn containing 
the number g(N, S) which gives the degeneracy of spin 
states. It is found 22 that 



g(N,s) 



N 



( N/2 - S 



N 

N/2 -5-1 



(14) 



Specifically for N — 4 particles, there is one spin eigen- 
function with S = 2, three with 5 = 1, and two with 
5 = 0. In general the spin part of the EXD wave func- 
tions involves a linear superposition over all the degener- 
ate spin eigenfunctions for a given 5. 

In the case of zero magnetic field and for a small num- 
ber of particles, one can find compact expressions that 
encompass all possible superpositions. For example, for 
N = 4 and 5 = 0, S z = one hasM 



mi) + (~ cos sin 0) | UU) 



X 00 = \/gSlEf 



-(-COS0 + 
-(-COS0 + 

+ (icos6> - 



— sin0)|itT4.> 



^sin0)|mT) 



(15) 



where the parameter 8 satisfies — ir/2 < 8 < tt/2 and is 
chosen such that 8 = corresponds to the spin function 



FIG. 2: Energy spectra (as a function of the magnetic field 
B) for N = 4 electrons in a double quantum dot with interdot 
separation d = 30 nm. Case of weak interelectron repulsion 
corresponding to GaAs (k = 12.5). The calculation were done 
for the case S z — and the Zeeman term was neglected. In 
this case all states with the same total spin S and different 
spin projections S z are degenerate. The effect of the Zeeman 
term can be easily added. Remaining parameters: e b = 0.5, 
huio = 5.1 meV, m* = 0.07m e . For all figures, the parameters 
hiVo and m* are kept the same (see Section [TTJ) . Energies are 
referenced to Nh^/uig + wf/4, where lu c = eB/(m*c) is the 
cyclotron frequency. 
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K = 6 
d=30 nm 




P=+1 ,S=0 



p=+1,S=2 
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FIG. 3: Energy spectra (as a function of the magnetic field B) 
for N = 4 electrons in a double quantum dot with interdot 
separation d = 30 nm. Case of intermediate interelectron 
repulsion (« = 6). The interdot barrier corresponds to e b = 
0.5. Energies are referenced to NTi^Jljq + oj|/4, where lj c = 
eB/(rn*c) is the cyclotron frequency. 



with intermediate two-electron spin S±2 = and three- 
electron spin 5123 = 1/2; whereas 8 = ±ir/2 corresponds 
to the one with intermediate spins S12 = 1 and ^123 = 
1/2. 
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For N = 4 and 5 = 1, 5, = one has: 



<TlO 



sm sin </3 



sin t) sm <p — 
sin 9 cos 
sin cos if 
sin sin ip + 
sin sin ip + 



^-sin0cos<£- icosWj | j j ] ) 
— sin 9 cosip + -cosV) \\\\\) 



^sinflsin^- ^costf) UTIT) 



sin6»cos(y3) Ittll) 
sin0cos^)|||TT>, (16) 



where the parameters 9 and tp satisfy — ir/2 < 9 < tt/2 
and — 7r/2 < ip < ir/2. Three independent spin functions 
with definite intermediate two-electron, 5i2, and three- 
electron, S*i23, spin values correspond to the 9 and ip 
values as follows: for S12 = and 6*123 = 1/2, 9 = 0; for 
S12 = 1 and S123 = 1/2, 9 = ±ir/2 and <p = 0; and for 
S12 = 1 and 5123 = 3/2, 9 = ±ir/2 and p = ±ir/2. 

Finally, for N = 4 and S = 2, S g = (maximum 
polarization) case, one has: 

<%20 = 

inn) + inn) + um> + inn) + inn) + inn) 



V6 



(17) 



V. RESULTS: ENERGY SPECTRA 



The excitation spectra as a function of the applied 
magnetic field for four electrons in a double QD with 
interdot distance d = 2x2 = —2x± = 30 nm and no volt- 
age bias between the dots [hi = hi = 0, see Eq. (TJJ)] are 
plotted for three different values of the interelectron re- 
pulsion strength, i.e., weak [k = 12.5 (GaAs); see Fig. [5], 
intermediate (k = 6; see Fig. [3]), and strong (k = 2; see 
Fig. 0]) Coulomb repulsion. The interdot barrier param- 
eter was taken as e b = 0.5 (because hi = /12 = 0, one has 
e i =e 2 = eb 'i see Section [TT1 for the definitions). In all 
cases, we calculated the eight lowest energy levels. 

We observe that the lowest six levels form a band that 
separates from the rest of the spectrum through the open- 
ing of a gap. This happens already at a relatively weak 
interelectron repulsion, and it is well developed for the 
intermediate case (k = 6). It is of interest to note that 
the number of levels in the band (six) coincides with the 
total number of spin eigcnfunctions for N — 4 fermions, 
as can be seen from the branching diagram displaying 
the spin degeneracies. In particular, there is one level 
with total spin S = 2 (and parity V xy — l)i three levels 
with total spin 5=1 (two with V xy = 1 and one with 
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FIG. 4: Energy spectra (as a function of the magnetic field 
B) for N = 4 electrons in a double quantum dot with interdot 
separation d — 30 nm. Case of strong interelectron repulsion 
(k = 2). The interdot barrier corresponds to e b = 0.5. Ener- 
gies are referenced to NH^/lOq + wf/4, where u c = eB/(m*c) 
is the cyclotron frequency. 
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FIG. 5: Energy spectra (as a function of the magnetic field 
B) for N = 4 electrons in a double quantum dot with in- 
terdot separation d — 60 nm. Case of weak interelectron 
repulsion corresponding to GaAs (k = 12.5). The interdot 
barrier corresponds to e b = 0.5. Energies are referenced to 
NH^luq + oj^/4, where uj c = eB/(m*c) is the cyclotron fre- 
quency. 



V xy = —1), and two levels with total spin 5 = (one 



with V X y = 1 and the second with V xy = — 1). All these 
six levels approximately cross at one point situated at 
about Bk3.5T for k = 12.5 and B » 2.2 T for k = 6. 

The trends associated with the opening of a gap and 
the formation of a six-state low band appear further rein- 
forced for the larger interdot distance of d = 60 nm (dis- 
played in Figs. [S] — [7] for the three values of the dielectric 
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constant re — 12.5, 6, and 2, respectively). It is remark- 
able that the six lower curves cross now at a sharply 
defined point (situated at B sa 3.3 T for re = 12.5 and 
B ps 2.1 T for re = 6. The six curves demonstrate ad- 
ditional near degeneracies regrouping approximately to 
three curves before and after the crossing point, which 
results in a remarkable simplification of the spectrum. 



For strong repulsion (re — 2), all six states in the 
low band are practically degenerate for both distances 
(d = 30 nm ; see Fig. 2] and d = 60 nm; see Fig. [7]). 
This is a consequence of the formation of a near-rigid 
Wigner molecule (WM) with strongly localized electrons. 
Namely, the overlap between the orbitals of localized elec- 
trons are practically zero (see, e.g., Ref. [23h . yielding a 
vanishing exchange, 38 and thus all six possible spin mul- 
tiplicities become degenerate in energy. Furthermore the 
physical picture of a near-rigid Wigner molecule suggests 
that the energy gap to the next band of states corre- 
sponds to excitation of the lowest stretching vibrational 
mode of the 4-electron molecule. 



It is natural to anticipate at this point that the above 
behavior at low B can be generalized to an arbitrary 
number of electrons TV in a double QD. Namely, as the 
strength of the interelectron interaction increases, a low- 
energy band comprising all possible spin multiplicities 
will form and it will become progressively well separated 
by an energy gap from the higher excitations. For ex- 
ample, for TV = 6, an inspection of the branching di- 
agram in Fig. [1] leads us to the prediction that there 
will be 20 states in this low-energy band. A similar be- 
havior emerges also in the case of a single, but strongly 
anisotropic quantum dot; indeed a low-energy band of 
three states (see the branching diagram in Fig. [T|) has 
been found for TV = 3 electrons in Ref. US 



It is of interest to contrast the above behavior of the 
excitation spectra in a double QD with that of an TV- 
electron circular dot. Specifically, in the circular QD, 
large inctclcctron repulsion leads to formation of a near- 
rigid rotating Wigner molecule that exhibits a rigid mo- 
ment of inertia. Then the states inside the low-energy 
band (two states for TV = 2, three for TV = 3, six for 
TV = 4, etc.) do not become degenerate in energy, but 
form an yrast rotational band 39 specified by L 2 J2Jq, 
where L is the total angular momentum and Jq is the 
classical moment of inertia. We note that the energy 
splittings among the yrast rotational states are much 
smaller than the vibrational energy gap in circular dots 
associated with the quantum of energy \/?>Tiujq of the 
stretching (often referred to as breathing) mode of the 
polygonal-ring configuration of the quasiclassical Wigner 
molecule j 40 i 41 ' 42 



VI. RESULTS: ELECTRON DENSITIES 

The electron density is the expectation value of the 
one-body operator 

N 

p(r) = (18) 

that is: 

p(r) = (^ X 9 D |p(r)|^ D ) 

= j2 c f c J^i\p(^j)- ( 19 ) 

Since p(r) is a one-body operator, it connects only 
Slater determinants and that differ at most by 
one spin orbital Xj ( r ) j f° r the corresponding Slater rules 
for calculating matrix elements between determinants for 
one-body operators in terms of spin orbitals, see Table 
2.3 in Ref. IIS 

In Figs. |HJa-f), we display (for the aforementioned 
three strengths of interelectron repulsion) the ground- 
state electron densities for for TV = 4 electrons in the 
case of a double dot at zero magnetic field with inter- 
dot separations d = 30 nm (left column) and d = 60 nm 
(right column). 

For the weak interaction case (re = 12.5) at B = 0, the 
electron densities do not exhibit clear signatures of for- 
mation of a Wigner molecule for either interdot distance, 
d = 30 nm [Fig. [SJa)] or d = 60 nm [Fig. Mb)]- The 
Wigner molecule is well formed, however, in the case of 
the intermediate Coulomb repulsion [re = 6; see Figs.EJc- 
d)]. One observes indeed four humps that correspond to 



K = 6 d=60 nm 




12 3 4 

B(T) 

FIG. 6: Energy spectra (as a function of the magnetic field B) 
for TV = 4 electrons in a double quantum dot with interdot 
separation d = 60 nm. Case of intermediate interelectron 
repulsion (re = 6). The interdot barrier corresponds to e b — 
0.5. Energies are referenced to NTi^/ujq + where lj c = 

eB/(m*c) is the cyclotron frequency. 



7 



97 



> 



95 





■ -».'.s=o ...» 


,3==* 

, : -.:-.-.--.*' ! '' 


'"■ ■- ..-**** 

K = 2 d=60 nm 






p=-i,s=i 


^^^=+1^=1 ^ 

" P=-1 ,s=l_ 



B(T) 



FIG. 7: Energy spectra (as a function of the magnetic field B) 
for N = 4 electrons in a double quantum dot with interdot 
separation d — 60 nm. Case of strong repulsion (k = 2). 
The interdot barri er correspo nds to e b = 0.5. Energies are 
referenced to Nh^Juj'^ + Uf/i, where u c = eB/(m*c) is the 
cyclotron frequency. 

the four localized electrons; they are located at (±34.88 
nm, ±13.13 nm) in the d = 60 nm case. In the case of 
strong Coulomb repulsion (k = 2) and for the same inter- 
dot distance d = 60 nm, the electrons are further local- 
ized as can be seen from Fig. [Hff ); the four humps occur 
now at (±39.86 nm, ±21.02 nm). The Wigner molecule 
is also well formed in the the strong-repulsion and d = 30 
nm case, as can be seen from Fig.|SJe), with the localized 
electrons located at (±29.28 nm, ±21.11 nm). 

VII. RESULTS: SPIN-RESOLVED 
CONDITIONAL PROBABILITY 
DISTRIBUTIONS AT B = 

A. Definitions 




d=30 nm d=60 nm 

FIG. 8: Electron density at B = for the ground state (with 
S = 0, S z = 0, parity P xy — 1) for N = 4 electrons in a 
double quantum dot with interdot separations d = 30 nm 
(left column) and d = 60 nm (right column). The top (a-b), 
middle (c-d), and bottom (e-f) rows correspond to ft = 12.5 
(weak repulsion), k — 6 (intermediate repulsion), and k = 2 
(strong repulsion), respectively. Ground-state energies: (a) 
E = 27.609 meV, (b) E = 21.572 meV, (c) E = 49.217 meV, 
(d) E = 39.799 meV, (e) E = 111.361 meV, (f) E = 94.516 
meV (compare Figs. [3] — EI) - The interdot barrier corresponds 
to e b = 0.5. Distances in nm. Vertical axis in arbitrary units 
(with the same scale for all six panels). 



In the regime corresponding to a well-defined Wigner 
molecule, the electron densities (see Sect. IVI[) are charac- 
terized by four humps that reflect the localization of the 
four electrons in the double quantum dot. Such charge 
densities do not provide any information concerning the 
spin structure of each EXD state. In fact, all six EXD 
states in the lower band exhibit very similar four-humped 
electron densities. 

The spin configurations associated with a given (5, S z ) 
EXD state in the WM regime can be explored with the 
help of the spin-resolved two-point anisotropic correla- 
tion function defined as: 

P CT(To (r,ro) = 

<$* XD | £ 6(r - r0*(r - r,-) W^, \*%?), (20) 



with the EXD many-body wave function given by equa- 
tion ©. 

Using a normalization constant 

W(cr,cr ,r ) = J P CTCTo (r,r )dr, (21) 

we further define a related conditional probability distri- 
bution (CPD) as 

V aao {v,r ) = P (7(To (r,r )/M(a,a ,r ), (22) 

having the property J "P CT(To (r, r )dr = 1. The spin- 
resolved CPD gives the spatial probability distribution 
of finding a second electron with spin projection a under 
the condition that another electron is located (fixed) at 
ro with spin projection <jq\ <J and cjq can be either up (f) 
or down Q). 
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FIG. 9: CPDs V n at B = for several EXD states with 

5 = 0, S z — 0, and parity P xy = 1 of N = 4 electrons in 
a double quantum dot with interdot separations d — 60 nm 
(a-c) and d = 30 nm (d). Case of strong Coulomb repulsion 
(ft = 2) with an interdot barrier e b = 0.5 (a,c-d) and e b = 

6 (b). Panels (a-b,d) correspond to ground states. Panel 
(c) corresponds to the excited second S = state for the 
same parameters as in panel (a) (see Fig. [7] and the branching 
diagram in Fig. [l}. Energies: (a) E = 94.516 meV, (b) E = 
96.811 meV, (c) E = 95.017 meV, and (d) E = 111.361 meV 
(compare Figs. [4]and[7]|. Distances in nm. Vertical axis in 
arbitrary units (with the same scale for all panels in Figs. 
[9] — lll|) . The fixed point is located at the maximum of the 
hump in the lower-left quadrant of the corresponding electron 
density, i.e., at ro =(—40 nm, —21 nm) for panels (a-c) and 
ro =( — 29 nm, —19 nm) for panel (d). 



Since T CT er ( r , r o) is a two-body operator, it connects 
only Slater determinants tyf and that differ at most 
by two spin orbitals Xji ( r ) and Xj 2 ( r )> f° r the correspond- 
ing Slater rules for calculating matrix elements between 
determinants for two-body operators in terms of spin or- 
bitals, see Table 2.4 in Ref. Km 



B. Examples of S = 0, S z = EXD states 

For each charge density corresponding to a given state 
of the system, one can plot four different spin-resolved 
CPDs, i.e., "Pjj, Vfi, Vii, and Vn- This can potentially 
lead to a very large number of time consuming compu- 
tations and an excessive number of plots. For studying 
the spin structure of the S = 0, S z = states at B = 0, 
however, we found that knowledge of a single CPD, taken 
here to be Vn (see Fig. [5]), is sufficient in the regime of 
Wigner-molecule formation. Indeed, the specific angle 
9 specifying the spin function Aoo Eq. (|15[) correspond- 
ing to the CPDs portrayed in Fig. [5] can be determined 
through the procedure described in the following: 

We designate with roman indices /, II, III, and TV 
the four quadrants of the (x, y) plane, starting with the 
upper left quadrant and going clockwise [see Fig. [DJa)]. 
In the case of a 4e Wigner-molecule, a single electron is 
localized within each quadrant. The same roman indices 
designate also the positions of the localized electrons in 
each of the six Slater determinants (e.g., | TTIDi I fiti)> 
etc.) that enter into the spin function Xqq in Eq. (|T5]) . We 
take always the fixed point to correspond to the fourth 
(IV) quadrant [bottom left in Fig. OJa)]. An inspection 
of Eq. ([15)) shows that only three Slater determinants 
in X ao contribute to V^, namely | TTU), I TIT!), and 
| ITTI); these are the only determinants in Eq. ([Pj]) with 
a down spin in the 4th quadrant. From these three Slater 
determinants, only the first and the second contribute 
to the conditional probability 11^(1) of finding another 
electron with spin-up in quadrant I: this corresponds to 
the volume under the hump of the EXD CPD in quadrant 
I [see, e.g., the hump in Fig. [SJa)]. Taking the squares 
of the coefficients of | TTII) and | TITI) m Eq. ([T5)l . one 
gets 



To calculate P aao (r, r ) in Eq. (|2"U)) . we use a sym- 
metrized operator 

f aao (r,r ) = 
^ [S(r - ri)S(r - r j )S aai 6 aorT .+ 



S(r - rj)S(r - r i )5 . (T ., <W*] , 



yielding 



(23) 



n n (/)oc 




(25) 



Similarly, one finds that only | TTII) an( i I ITTI) con- 
tribute to II||(/J), and that 



n n (jj)« 




(26) 



Integrating under the humps of the EXD CPD in 
quadrants / and //, we determine numerically the ratio 



= V Cf C q j(^f \f 1*^). (24) n u (/)/n u (/7), which allows us to specify the absolute 



/../ 



value of 9 (within the interval -90° < 9 < 90°) via the 
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expressions in Eqs. |25J and ([26]) . The restriction to the 
absolute value of 9 is a result of the squares of the sine 
and cosine entering in II ^ (/) and 11^(11). To obtain the 
actual sign of 9, additional information is needed: for ex- 
ample the ratio H^i(I)/H^(III) can be used in a similar 
way, where 



U n {III) tx (^cos( 



S=1 



\ cost 



-sin0 



sin ( 



(27) 



Using the method described above, we find that 9 « 
—60° for the EXD ground state at d = 60 nm (larger 
interdot distance) and k = 2 [strong repulsion; see Fig. 
OJa)], and the corresponding spin function simplifies to 

4o = -~i ttu>+|i uu>+~i uu)-~i urn- (28) 

Remarkably, increasing the interdot barrier from e b = 
0.5 [Fig. to e b = 6 [Fig. M^)}, while keeping the 

other parameters constant, does not influence much the 
composition of the associated spin function, which re- 
mains that given by Eq. (|28[) . This happens in spite 
of the visible change in the degree of localization in the 
electronic orbitals, with the higher interdot-barrier case 
exhibiting a sharper localization. 

In Fig. [9](c) , we display the Vn CPD for an excited 
state with S = 0, S z = (having P xy = 1 and energy 
E = 95.017 meV), with the remaining parameters being 
the same as in Fig. HJa). For this case, following an 
analysis as described above, we found the angle 9 « 30°, 
which is associated with a spin function of the form 



(2) 



1 
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Ull) 



i 



TITI) 



l 



Tilt) 



2V3 2V3 V3 

' ' iTT4) + -^l 4T4T) + t4| Utt). 



2V3 



2V3 



(29) 



and 



We note that the spin functions in Eqs. 
are orthogonal. 

In Fig. [5{d), we display the Pn CPD for the ground 
state with S = 0, S z — (having P xy = 1 and energy 
E = 111.361 meV) and for the shorter interdot distance 
d = 30 nm. For this case, we found an angle 9 ~ —63.08°, 
which corresponds to the following spin function: 



y(3) _ 



-0.5148| TTU) + 0.4838| TITI) + 0.031| TUT) 
0.031| |TU) + 0.4838| ITU) - 0.5148| ||H). 

(30) 



From a comparison of the above result with that for 
the larger d = 60 nm [see Eq. (f2"5|) ]. we conclude that the 
difference in interdot distance results in a slight variation 
of the spin functions. 



d=60 nm 




FIG. 10: CPDs at B = for excited EXD states with S = 
1,S Z = 0, and parity P xy = 1 of N = 4 electrons in a double 
quantum dot at the larger interdot separation d — 60 nm (a- 
b) and the shorter interdot separation d — 30 nm (c). Panels 
(a) and (c) display Vn CPDS (down- up), while panel (b) 
displays a different Vyy CPD (down-down), but for the same 
state as in (a) . Case of strong Coulomb repulsion (k = 2) with 
interdot barrier e ft = 0.5. Energies: (a-b) E — 94.757 meV, 
and (c) E — 111.438 meV (compare Figs. |3]and[7|. Distances 
in nm. Vertical axis in arbitrary units (with the same scale 
for all panels in Figs. [9] — lll[) . The fixed point is located 
at the maximum of the hump in the lower-left quadrant of 
the corresponding electron density, i.e., at ro =(—40 nm, —21 
nm) for panels (a-b) and ro =(—29 nm, —19 nm) for panel 
(c). Note that this is a case with 5 = 1; the previous Fig. [9] 
displayed S = cases. 



C. Examples of S = 1, S z = EXD states 

In this section, we turn our attention to partially po- 
larized EXD states with S = 1. 

In Fig. [TOta), we display the T n CPD at B = for 
an excited state with S = 1,S Z = 0, parity P xy = 1, 
and energy E — 94.757 meV, at the larger interdot 
separation d = 60 nm. Again we consider the case of 
strong Coulomb repulsion (k — 2) with an interdot bar- 
rier e b = 0.5. The corresponding spin function Aio [Eq. 
(fTo) ] depends on two different angles 9 and <j), and one 
needs at least two different CPDs for determining their 
specific values. For this purpose, we display also the Vn 
CPD for the same state in Fig. riWb) . 

The specific values of 9 and tfi associated with the 
CPDs in Figs. [TOT ai and [TOT bi can be determined 
through the ratios U n (I)/U n (II) and U n (I)/U n (III) 
[associated with Fig. [TPT a)] and Hii(I)/Ili[(II) and 
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n u (7)/II u (777) [associated with Fig. [T0t»]. where 



B=2T 



n n (7) cx -sin 2 6>shV>T — 
. V2 



5 1 
sin 2 9 cos 2 d> + - cos 2 i 



■ sin ft sin <z> cos <p H = sin ft cos ft sin ( 

6 V6 



12 
1 



sin ft cos ft cos ( 



(31) 



5 1 
n-f; (77) oc ^ sin 2 9 sin 2 0+^2 sin2 ^ cos2 ^ + 4 cos2 ' 



V2 



1 



■ sin ft sin cos = sm ft cos ft sin < 

6 V6 



sm ft cos ft cos 0, 



' ' 2 °~- 2 J- 1. i sin 2 ft cos 2 . 



(32) 



n Ti (777) cx -sin 2 ft sin 
o 



\[2 1 

— — - sin 2 6* sin (/) cos (/) + — cos 2 9, (33) 




— sm ft sm c 
6 



Ha (77) oc I \/ -sm ft sm.. - 



n u (777) cx l \/-smftsm, 



sin ft cos 6 — - cos t 
12 ^2 



(34) 



sm ft cos 




— sm ft cos < 



(36) 



Using Eqs. (T5T|) — (|3"6"1) and the numerical val- 
ues of the ratios H. n (7) /Tl n (77) and II n (7)/II n (777) 
n U (7)/n u (77) and n u (7)/n u (777) (specified via a 
volume integration under the humps of the EXD CPDs), 
we determined that 9 = —45° and sin0 = — -v/2/3, 
cos0 = (i.e., (f) ps -54.736°). Thus, the corre- 

sponding spin function reduces to the simple form 

^io = ^|itm>-^/|imt>. (37) 

In Fig. [TDtc), we display the 7> u CPD at B = for a 
similar excited state as in Fig. nW a) (with S = 1,S Z = 0, 
parity P XJ/ = 1, and energy 75 = 111.438 meV) of N = 4 
electrons at the shorter interdot separation d — 30 nm. 
Here too we consider the case of strong Coulomb repul- 
sion (k = 2) with interdot barrier e b — 0.5. We note that 
the localization of electrons is stronger for the larger in- 
terdot distance [compare Fig.fTOTai with Fig. flQf c)]. This 
difference, however, does not influence the coefficients en- 
tering into the associated spin function, which we found 
to remain very close to the specific form in Eq. (|37[) . 




lower energy 



higher energy 



FIG. 11: V n CPDs for the EXD states at B = 2 T with S = 0, 
S z ~ 0, and parity P xy = 1 of N = 4 electrons in a double 
quantum dot at the larger interdot separation d — 60 nm. 
(a) The lower energy of the two S — states (see branching 
diagram in Fig. [1]). (b) Higher energy S — state. Case 
of strong Coulomb repulsion (k — 2) with interdot barrier 
e b = 0.5. Energies: (a) E = 94.605 meV and (b) E = 95.047 
meV (compare Fig. [7]). Distances in nm. Vertical axis in 
arbitrary units (with the same scale for all panels in Figs. [5] 
— 11 ip . The fixed point is located at ro =(—40 nm, —21 nm). 



VIII. RESULTS: SPIN-RESOLVED 
CONDITIONAL PROBABILITY 
DISTRIBUTIONS 

In Fig. [TTJ we display EXD CPDs at a finite value of 
the magnetic field, and precisely at B = 2 T, for the 
two states of the low-energy band with S — 0, S z = 
(at the larger interdot separation d = 60 nm and strong 
interelectron repulsion k = 2). This value of B was cho- 
sen to lie beyond the crossing point for the six states of 
the low-energy band (which happens at B ~ 1 T; see 
Fig. [7]). Comparison with the CPDs of the correspond- 
ing states at zero magnetic field [sec Figs.[5Ja) and^c)] 
shows that the spin structure of the associated Wigner 
molecule varies rather slowly with the increasing mag- 
netic field in the range < B < 2.5 T. 

Following the height of the humps in the left upper 
quadrants, one observes that the CPD in Fig. lTlT a) (case 
of lower- energy state at B — 2 T with S = and P xy = 1) 
corresponds to that of Fig. [5{a) (case of lower-energy 
state at B = with 5 = and P xy = 1). Similarly, 
the CPD in Fig. rTJTb) at 7? = 2 T {higher- energy state) 
corresponds to that of Fig. M,c) at B — (higher- energy 
state). From these results, we concludes that the two 



states with S = and P xy — 1 do not really cross at 
the 'crossing' point at B ~ 1 T. In reality, this point is 
an anticrossing point for these two states, although the 
anticrossing gap is too small to be seen with the naked 
eye. This behavior agrees with that expected from states 
having the same quantum numbers. We checked that a 
similar observation applies for the two other states in the 
low-energy band having the same quantum numbers, i.e., 
those having S = 1 and P xy = — 1. 
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IX. DISCUSSION 
A. Finite Heisenberg spin clusters 

In IVII| using the spin-resolved CPDs, we showed that 
the EXD many-body wave functions in the Wigner- 
molecule regime can be expressed as a linear superpo- 
sition of a small number of Slater determinants and that 
this superposition exhibits the structure expected from 
the theory of many-body spin functions. This finding 
naturally suggests a strong analogy with the field of nano- 
magnets and quantum magnetism, usually studied via 
the explicitly spin-dependent model effective Hamilto- 
nian known as the Heisenberg Hamiltonian , 12 ' 13 i 43 given 
by: 

T~£'h — ^ JijSjSj — Q Si, (38) 

i,3 » 

where Jy are the exchange integrals between spins on 
sites i and j. Even in its more familiar, simplest form 

n H = Jj2 s *- s 3> ( 39 ) 

that is that of the spin- 1/2 Heisenberg antiferromagnet 
with nearest-neighbor interactions only, it is well known 
that the zero-temperature (at B = 0) solutions of Hamil- 
tonian ([55]) involve radically different forms as a function 
of the geometry, dimensionality, and size. 

It is natural to compare the EXD spin functions de- 
termined in Section IVIII with well known solutions of the 
Heisenberg Hamiltonian [Eq. ([55]) ] when the four spins 
are located on four sites arranged in a perfect square i 13 i 14 
(The perfect-square arrangement arises 4 ^ also in the case 
of formation of a four-electron Wigner molecule in a sin- 
gle circular quantum dot.) In this case, the ground state 
of H.h is the celebrated resonating valence bond (RVB) 
stat o 13 ' 14 which forms the basic block in many theoret- 
ical approaches aiming at describing high-temperature 
superconductors^ The RVB state has quantum num- 
bers S = 0, S z = and is given b y 13 i 14 

"*°° VB = 273 lTTU>+ 2vl lTUT>+ 273 IUTT> 

+ 27f um> ~ 7! 1 nn> ~ 7! lmT> ' 

(40) 

(2) 

Although the excited-state EXD Xq Q in the quantum- 
double-dot case portrayed in Fig.[9lc) appears to be sim- 
ilar to the RVB X^ B [Eq. (3DJ)], they are not equal. 
Indeed the coefficients of the pair of Slater determimants 
| tltl) an£ l I itlf) have been interchanged with those of 

| |TT4> and | TUT). 

Similar observations apply also to the 5 = and S z 
( j remaining states that are orthogonal to Xqq [see X^ 



in Eq. (f28|) : case of double quantum dot] and to A 1 ^ 6 
(case of a perfect square). The latter is given byi 3 - 

x „ exci = _i ] TTU)+ i| TUT)+ I| im)~\ un>- 

(41) 

In particular, one finds 

y(2) _ L vRVB are,cxci f AO\ 

^oo - _ 2 00 2~ 00 

and 

v (l) _ \/3 v p{VB | 1 v squarc,cxci ( A1\ 

<*oo — ^—^vq + 2 00 ' ^ ' 

We note that the differences in the Abo spin functions 
between the DQD case (corresponding to a parallelo- 
gram) and the perfect-square case are also reflected in the 
V n CPDs. Indeed the CPDs of the DQD exhibit equal- 
height humps along the smaller side of the parallelogram 
while those of the perfect-square configuration (and/or 
circular quantum dot) exhibit equal-height humps along 
a diagonal.— 

Neel antiferromagnetic ordering, where the average 
spin per site < S* >= ( — l) J+1 /2, is an important mag- 
netic phenomenon in the thermodymanic limit 1 - 4 " associ- 
ated with breaking of the total-spin symmetry. The finite 
size magnetic clusters discussed here exhibit a sharply 
different behavior in this respect. Indeed, as discussed 
in Ref. [3, the four-site Neel state is the single Slater 
determinant | HIT) ( or T-it-i))- It is clear that the total- 
spin conserving EXD functions Xqq are multideterminen- 
tal and have an average spin per localized electron (per 
site) < >= 0. 

We concur with Ref. [U that the phenomenon of Neel 
antifferomagnetism is not applicable to assemblies of few 
electrons. In the next section, we argue that the appro- 
priate concept for WM states is that of spin entangle- 
ment. 



B. Spin entanglement 

In the previous sections, we showed that the EXD wave 
functions in the regime of Wigner- molecule formation can 
be approximated as a superposition of a small number 
of Slater determinants corresponding to well structured 
spin functions; see, e.g., X$ in Eq. (28]). This is a 
great simplification compared to the initial EXD super- 
position [Eq. |7|)]. where the counting index is usually 
I > 500, 000. This reduction of the molecular EXD solu- 
tions to their equivalent spin functions (described in Sec- 
tion IVIip enables one to investigate their properties re- 
garding fundamental quantum behavior associated with 
quantum correlations and fluctuations beyond the mean 
field. 

The smallest number of Slater determinants contribut- 
ing to the spin functions X's with S z = is two. We 
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note that for S z — a single determinant of four lo- 
calized spin orbitals cannot conserve the total spin, and 
thus mean-field approaches like DFT (which are associ- 
ated with a single determinant) are unable to describe 
quantum fluctuations and entanglement (see below). In 
particular, we note that the spin functions ?Cs cannot be 
further reduced to simpler superpositions containing a 
smaller number of Slater determinants. As a result, they 
faithfully represent the extent of probabilistic quantum 
interconnection between the individual electrons partic- 
ipating in the system as described by an EXD solution. 
This quantum interconnection is widely referred to as 
entanglemen t 46 'i 47 i 48 and generates correlations between 
the physical observables of the individual electrons (e.g., 
spins), even though the individual electrons are spatially 
separated. This existence of quantum interconnectivity 
independently of spatial separation has made entangle- 
ment the central instrument for the development of the 
fields of quantum information, such as quantum telepor- 
tation, quantum cryptography, and quantum computing. 
With their localized electrons, Wigner molecules in quan- 
tum dots offer another physical solid-state nanosystem 
where entanglement may be realized and studied. 

The mathematical theory of entanglement is still devel- 
oping and includes several directions. One way to study 
entanglement is through the use of properly defined mea- 
sures of entanglement, e.g., the von Neumann entropy 
which utilizes the single-particle density matrix. Another 
way is to catalog and specify classes of entangled states 
that share common properties regarding multipartite en- 
tanglement. A well known class of A-qubit entangled 
states are the Dicke states ; 24 ' 25 ! 26 ' 27 which most often 
are taken to have the symmetric form: 

XBt e =i k ) (|n^_100...0)+Perm). (44) 

k 

Each qubit is a linear superposition of two single-particle 
states denoted by or 1, and the symbol 'Perm' stands 
for all remaining permutations. The or 1 do not have 
to be necessarily up or down 1/2-spin states. Two-level 
atoms in linear ultracold traps have already been used 
as an implementation of a qubit. Dicke states appear 
in many physical processes like superradiance and super- 
fluorescence. They can also be realized with photons, 
where the qubits correspond to the polarization degree 
of freedom!^ 7 . 

In the 1/2-spin case of fermions (e.g., for electrons), the 
Dicke states of Eq. (|4"4"|) correspond to a fully symmetric 
flip of k out of N localized spins. It is apparent that the 
four-qubit fully polarized (S — 2 with spin projection 
S z =0) EXD solution is reproduced by A20 01 Eq- (fT7| . 
and thus it is of the symmetric Dicke form (with k = 2) 
displayed above in Eq. (|4~4"|) . On the other part, the DQD 
EXD states (with S z = 0) studied in Section IVIII with 
S = and/or S = 1 represent a natural generalization 
of Eq. (|4"4"1) to the class of asymmetric Dicke states. We 
hope that our results will motivate experimental research 



aiming at the realization and control of such states in 
DQD electronic devices. 

Before leaving this section, we note that Dicke states 
with a single flip (k — 1) are known as W statesj 49 ' 50 
For N = 4 electrons, the latter states are related to EXD 
solutions with S z — ±1. For the connection between W 
states and EXD states for N = 3 electrons in anisotropic 
quantum dots, see Ref. [H. W states have already been 
realized experimentally using two-level ultracold ions in 
linear traps^ 



X. SUMMARY 

Extensive investigations of lateral double quantum 
dots containing four electrons were performed using the 
exact-diagonalization method, as a function of interdot 
separation, applied magnetic field, and strength of inter- 
electron repulsion. Novel quantum behavior was discov- 
ered compared to circular QDs, concerning both energy 
spectra and quantum entanglement aspects. Thus it is 
hoped that the present work will motivate further exper- 
imental studies on lateral DQDs beyond the two-electron 
casej 2 ^ 5 - 

Specifically it was found that, as a function of the mag- 
netic field, the energy spectra exhibit a low-energy band 
consisting of a group of six states, and that this number 
six is not accidental, but a consequence of the conserva- 
tion of the total spin and of the ensuing spin degenera- 
cies and supermultiplicities expressed in the branching 
diagram. These six states appear to cross at a single 
value of the magnetic field, and the crossing point gets 
sharper for larger interdot distances. As the strength of 
the Coulomb repulsion increases, the six states tend to 
become degenerate and a well defined energy gap sepa- 
rates them from the higher-in-energy excited states. 

The formation of the low-energy band is a consequence 
of the localization of the four electrons within each dot 
(with two electrons on each dot). The result is forma- 
tion (with increasing strength of the Coulomb repulsion) 
of a Wigner supermolecule, with the four localized elec- 
trons at the corners of a rectangular parallelogram. Using 
the spin-resolved pair-correlation functions, it was shown 
that one can map the EXD many-body wave functions to 
the spin functions associated with the four localized elec- 
trons. This mapping led us naturally to studying analo- 
gies with finite systems described by model Heisenberg 
Hamiltonians (referred to often as finite Heisenberg clus- 
ters) . It was found that the determination of the equiva- 
lent spin functions enables investigations concerning the 
entanglement properties of the EXD solutions. In partic- 
ular, it was shown that the formation of Wigner super- 
molecules generates strongly entangled states known in 
the literature of quantum information as A-qubit Dicke 
states^^ 2 ^! 
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FIG. 12: Single-particle spectra of a double quantum at B — 
plotted versus the distance d between two (identical) cou- 
pled QDs with a TCO confinement Ttw x i = fitj X 2 ~ huj y = 3 
meV and hi — hi — [see Eq. (JTJ)]. For all d's the bar- 
rier control parameters were taken as e\ = e b 2 = 0.5, i.e., 
the barrier height (depicted by the dashed line) varies as 
14(d) = Vo(d)/2. Molecular orbitals correlating the united 
(Vb = 0) and separated-dots limits are denoted along with the 
corresponding (on the right) single-QD states. Wave function 
cuts at y = along the a;-axis at several distances d (see ar- 
rows) corresponding to the lowest bonding and antibonding 
eigenvalues (solid and dashed lines, respectively) are displayed 
at the top. Energies in meV and distances in nm. 



Acknowledgments 

This work was supported by the US D.O.E. (Grant No. 
FG05-86ER45234). 



APPENDIX: SINGLE-PARTICLE STATES OF 
THE TWO-CENTER OSCILLATOR 

In this Appendix, we discuss briefly the energy spec- 
tra associated with the single-particle states of the two- 
center oscillator Hamiltonian given by Eq . JTJ) . We follow 
here the notation presented first in Ref. [32|. For further 
details, see Ref. (33- 

The calculated two-center oscillator single-particle 
spectrum for a double quantum dot made of two tunnel- 
coupled identical QDs (with fuv y = Tko x i = Tko X 2 = 3 
meV) plotted versus the distance, d, between the centers 
of the two dots, is given in Fig. [T3J In these calcula- 
tions, the height of the barrier between the dots varies as 
a function of d, thus simulating reduced tunnel coupling 
between them as they are separated; we take the barrier 
control parameter as e\ = e h 2 = 0.5. In the calculations 
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FIG. 13: Single-particle spectrum of the d = 70 nm (Tuj x i = 
Tijjj X 2 = hu) y — 3 meV, Vb = 2.43 meV, hi = hi = 0) double 
quantum dot versus B (in T). The hui c /2 (Ml = 0, lower line) 
and 3hui c /2 (Ml = 1, upper line) first and second Landau 
levels are given by the dashed lines. 



in this Appendix, we used GaAs values, m* = 0.067m e 
and a dielectric constant n = 12.9. For the separated 
single QDs (large d) and the unified QD (d = 0) limits, 
the spectra are the same, corresponding to that of a 2D 
harmonic oscillator (being doubly degenerate for the sep- 
arated single QDs) with a level degeneracy of 1, 2, 3, ... . 
In analogy with real molecules, the single-particle states 
in the intermediate region (d > 0) may be interpreted as 
molecular orbitals (MOs) made of linear superpositions 
of the states of the two dots comprising the DQ This 
qualitative description is intuitively appealing, though it 
is more appropriate for the weaker coupling regime (large 
d); nevertheless we continue to use it for the whole range 
of tunnel-coupling strengths between the dots, including 
the strong coupling regime where reference to the states 
of the individual dots is only approximate. Thus, for 
example, as the two dots approach each other, the low- 
est levels (n x , n y ) with n x = n y = on the two dots 
may combine symmetrically ( "bonding" ) or antisymmet- 
rically ("antibonding") to form [0,0;0] and [0,0;1] MOs, 
with the third index denoting the total number of nodes 
of the MO along the interdot axis (x), that is, 2n x +T, 
T = or 1; for symmetric combinations (X — 0), this 
index is even and for antisymmetric ones (X = 1), it is 
odd. Between the separated-single-QDs and the unified- 
QD limits, the degeneracies of the individual dots' states 
are lifted, and in correlating these two limits the num- 
ber of x-nodes is conserved; for example the [0,0;1] MO 
converts in the unified-QD limit into the (1,0) state of a 
single QD, the [1,0;2] MO into the (2,0) state, and the 
[0,1;1] MO into the (1,1) state (see Fig.HU). Note that 
MOs of different symmetries may cross, while they do 
not if they are of the same symmetry. 

In a magnetic field, the TCO model consitutes a gener- 
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alization of the Darwin-Fock models for non-interacting 
electrons in a single circular QD. The single-particle spec- 
tra for the DQD (d = 70 nm, Vf, = 2.43 meV) in a mag- 
netic field (B) are shown in Fig. 02] (here we neglect the 
Zeeman interaction which is small for our range of B 
values with g* = —0.44 for GaAs). The main features 
are: (i) the multiple crossings (and avoided crossings) as 
B increases, (ii) the decrease of the energy gap between 
levels, occurring in pairs (such as the lowest bonding- 
antibonding pair), portraying an effective reduced tun- 



nel coupling between the QDs comprising the DQD as B 
increases, (iii) the "condensation" of the spectrum into 
the sequence of Landau levels (Ml + l/2)huj c , Ml = 0, 
1, 2, ... (the Ml — and Ml = 1 bands are depicted, 
respectively, by the lower and upper dashed lines in Fig. 
n"3|) . This is similar to the behavior of the single-particle 
Darwin-Fock spectrum for harmonically confined elec- 
trons in a circular QD^ (note however that the geometry 
of the DQD is non-circular and deviates from a simple 
harmonic confinement). 



* Electronic address: Constantine.Yannouleas@physics.gatech.edu 
t Electronic address: Uzi.Landman@physics.gatech.edu 

1 L.P. Kouwenhoven, D.G. Austing, and S. Tarucha, Rep. 25 
Prog. Phys. 64, 701 (2001). 

2 R. Hanson, L.P. Kouwenhoven, J.R. Petta, S. Tarucha, and 26 
L.M.K. Vandersypen, Rev. Mod. Phys. 79, 1217 (2007). 

3 P.A. Maksym, H. Imamura, G.P. Mallon, and H. Aoki, J. 27 
Phys.: Conden. Matter 12, R299 (2000). 

4 S.M. Reimann and M. Manninen, Rev. Mod. Phys. 74, 28 
1283 (2002). 

5 C. Yannouleas and U. Landman, Rep. Prog. Phys. 70, 2067 29 
(2007). 

6 D. Loss and D.P DiVincenzo, Phys. Rev. A 57, 120 (1998). 

7 G. Burkard, D. Loss, and D.P. DiVincenzo, Phys. Rev. B 30 
59, 2070 (1999). 

8 R. Woodworth, A. Mizel, and D.A. Lidar, J. Phys.: Con- 
dens. Mat. 18, S721 (2006). 31 

9 S. Tarucha, D.G. Austing, T. Honda, R.J. van der Hage, 

and L.P. Kouwenhoven, Phys. Rev. Lett. 77, 3613 (1996). 32 

10 M. Ciorga, A.S. Sachrajda, P. Hawrylak, C. Gould, P. Za- 
wadzki, S. Jullian, Y. Feng, and Z. Wasilewski, Phys. Rev. 

B 61, R16315 (2000). 33 

11 M. Avinun-Kalish, M. Heiblum, O. Zarchin, D. Mahalu, 
and V. Umansky, Nature 436, 529 (2005). 

12 P.V. Hendriksen, S. Linderoth, and P.A. Lindgard, Phys. 
Rev. B 48, 7259 (1993). 

13 S. Haas, in Lectures on the Physics of Strongly Correlated 
Systems XII, edited by A. Avella and F. Mancini (AIP 
Conf. Proceedings Vol. 1014, Melville, New York, 2008). 

14 P. Fazekas, Lecture notes on electron correlation and mag- 34 
netism (World Scientific, Singapore, 1999) 

15 M.M Tavernier, E. Anisimovas, F.M. Peeters, B. Szafran, 
J. Adamowski, and S. Bednarek, Phys. Rev. B 68, 205305 
(2003). 35 

16 M. Rontani, C. Cavazzoni, D. Bellucci, and G. Goldoni, J. 
Chem. Phys. 124, 124102 (2006) 

17 A. Ghosal, A.D. Giiclii, C.J. Umrigar, D. Ullmo, and H.U. 
Baranger, Phys. Rev. B 76, 085341 (2007). 36 

18 M. Helle, A. Harju, and R.M. Nieminen, Phys. Rev. B 72, 
205329 (2005). ' 37 

19 L.X. Zhang, D.V. Melnikov, and J. P. Leburton, Phys. Rev. 
B 74, 205306 (2006). 

20 B. Szafran, and F.M. Peeters, Phys. Rev. B 71, 245314 
(2005). 

21 T. Chwiej and B. Szafran, Phys. Rev. B 78, 245306 (2008). 

22 R. Pauncz, The Construction of Spin Eigenf unctions: An 
Exercise Book (Kluwer Academic/Plenum Publishers, New 38 
York, 2000). 

23 C. Yannouleas and U. Landman, Phys. Rev. Lett. 82, 5325 39 



(1999); ibid. 85, E2220 (2000). 

R.H. Dicke, Phys. Rev. 93, 99 (1954). 

F. Verstraete, J. Dehaene, B. De Moor, and H. Verschelde, 

Phys. Rev. A 65, 052112 (2002). 

J.K. Stockton, J.M. Geremia, A. C. Doherty, and H. 
Mabuchi, Phys. Rev. A 67, 022112 (2003). 
J. K. Korbicz, O. Giihne, M. Lewenstein, H. Haffner, C. 
F. Roos, and R. Blatt, Phys. Rev. A 74, 052319 (2006). 
Yuesong Li, C. Yannouleas, and U. Landman, Phys. Rev. 
B 76, 245310 (2007). 

A three-dimensional magnetic-field-free version of the 
TCO has been used in the description of fission in metal 
clusters 3 "^ and atomic nuclei.— 

C. Yannouleas and U. Landman, J. Phys. Chem. 99, 14577 
(1995). C. Yannouleas, R.N. Barnett, and U. Landman, 
Comments At. Mol. Phys. 31, 445 (1995). 
J. Maruhn and W. Greiner, Z. Phys. 251, 431 (1972). C.Y. 
Wong, Phys. Lett. 30B, 61 (1969). 

R.N. Barnett, C.L. Cleveland, H. Hakkinen, W.D. 
Luedtke, C. Yannouleas, and U. Landman, Eur. Phys. J. 
D 9, 95 (1999). 

The two-dimensional TCO was used earlier to approximate 
the solutions of the many-body Hamiltonian in Eq. Q 
at the level of unrestricted Hartree-Fock (UHF) exhibit- 
ing breaking of total-spin symmetry, see Ref. [23|; C. Yan- 
nouleas and U. Landman, Eur. Phys. J. D 16, 373 (2001); 
Int. J. Quantum Chem. 90, 699 (2002); J. Phys.: Condens. 
Matter 14, L591 (2002). For a review of the necessary sub- 
sequent step of symmetry restoration, see Ref. 0. 
Due to the very long CPU times required for calculating 
surfaces, several CPDs at B = were calculated using K — 
32. We checked that this does not influence our analysis of 
spin functions compared to K = 50. 

R.B. Lehoucq, D.C. Sorensen, and C. Yang, ARPACK 
Users' Guide: Solution of Large-Scale Eigenvalue Prob- 
lems with Implicitly Restarted Arnoldi Methods (SIAM, 
Philadelphia, 1998). 

A. Szabo and N.S. Ostlund, Modern Quantum Chemistry 
(McGraw-Hill, New York, 1989). 

For the expression Xoo in Eq. (| 15ft . see Ch. 6.6 in Y. 
Suzuki and K. Varga, Stochastic Variational Approach 
to Quantum-Mechanical Few-Body Problems (Springer, 
Berlin, 1998). For the remaining two expressions for S z = 
0, i.e., Xio and X20 [Eqs. (Tl6jl and ([17)1 . respectively; see 
Ying Li, Ph.D. Dissertation, Georgia Institute of Technol- 
ogy (2009)]. 

C. Yannouleas and U. Landman, J. Phys.: Condens. Mat- 
ter 14, L591 (2002), in particular Fig. 3 therein. 
C. Yannouleas and U. Landman, Phys. Rev. B 69, 113306 



15 





(2004). 


48 


L. Amico, R. Fazio, A. Osterloh, and V. Vedral, Rev. Mod. 


40 


C. Yannouleas and U. Landman, Phys. Rev. Lett. 85, 1726 




Phys. 80, 517 (2008). 




(2000). 


49 


W. Diir, G. Vidal, and J.I. Cirac, Phys. Rev. A 62, 062314 


41 


V.A. Schweigert and F.M. Peeters, Phys. Rev. B 51, 7700 




(2000). 




(1995). 


50 


V. Coffman, J. Kundu, and W.K. Wootters, Phys. Rev. A 


42 


M.R. Geller and G. Vignale, Phys. Rev. B 53, 6979 (1996). 




61, 052306 (2000). 


43 


E. Manousakis, Rev. Mod. Phys. 63, 1 (1991). 


51 


H. Hiiffner, W. Hansel, C.F. Roos, J. Bcnhelm, D. Chek-al- 


44 


C. Shi, G.S. Jeon, and J.K. Jain, Phys. Rev. B 75, 165302 




kar, M. Chwalla, T. Korber, U.D. Rapol, M. Riebe, P.O. 




(2007). 




Schmidt, C. Becher, O. Giihne, W. Diir, and R. Blatt, 


45 


Ph. Anderson, Physics Today (April 2008) p. 8 


52 


Nature 438, 643 (2005). 


46 


E. Schrodinger, Proc. Cambridge Phil. Soc. 31, 555 (1935). 


V. Fock, Z. Phys. 47, 446 (1928); C.G. Darwin, Proc. Cam- 


47 


A. Einstein, B. Podolsky, and N. Rosen, Phys. Rev. 47, 




bridge Philos. Soc. 27, 86 (1930). 




0777 (1935). 







